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Abstract 


The  problem  of  diffraction  of  scalar  waves  by  an  infinite  conducting 
plane  with  a  slit  is  investigated.  Approximate  expressions  for  the  near 
and  far  fields,  taking  into  account  the  interaction  between  the  edges,  are 
derived  in  terms  of  the  well-known  solutions  for  the  field  produced  when 
an  isolated  conducting  half -plane  is  excited  by  a)  a  plane  wave,  and  b)  a 
line  source.  Results  of  numerical  calculation  are  given  for  the  case  of 
a  plane  wave  normally  incident  on  the  slit.  Twelve  values  of  slit  width 
reuiging  from  .96  to  2.5  wavelengths  are  considered.  For  each  slit  width 
the  far  diffracted  field  is  ccsnputed  at  ten-degree  intervals  in  our  approx- 
imation, and  is  compared  with  the  exact  solution  obtained  by  Mathieu  func- 
tions, and  with  the  results  obtained  when  interaction  is  neglected.  A 
comparison  of  transmission  coefficients  is  also  given.  It  is  found  that 
the  new  approximate  solution  agrees  well  with  the  exact  solution,  and  pro- 
vides a  significant  correction  to  the  non-interaction  solution.  The  accuracy 
increases  with  the  slit  width,  so  that  the  result  is  useful  in  the  range  where 
interaction  cannot  well  be  neglected  but  where  the  exact  solution  converges 
so  slowly  that  computation  is  impracticable,  A  brief  discussion  of  the  case 
of  line-source  excitation  is  included. 
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I,   Introduction 

Since  there  already  exist  numerous  investigations  of  the  diffraction 
of  scalar  waves  by  a  slit,  the  addition  of  another  item  to  the  literature  requires 
some  explanation.  The  purpose  of  the  present  note  is  to  give  an  elementary  and 
intuitive  account  of  the  diffraction  by  an  infinite  conducting  plane  with  a  slit. 
We  consider  specifically  the  case  of  electromagnetic  waves  whose  electric  field 
possesses  only  a  single  component,  namely  a  component  in  the  direction  parallel 
to  the  edges  of  the  slit.  The  results  are  also  applicable  to  the  acoustic  pro- 
blem, for  soft  screens.  It  is  found  that  the  results  for  large  slit  widths  can 
be  expressed  approximately  in  terms  of  simple  well-known  functions  and  that  these 
approximations  are  more  accurate  than  those  conveniently  obtainable  by  previous 
procedures. 

The  procedures  that  have  hitherto  been  proposed  and  applied  to  this 

problem  fall  into  three  categories.  First,  there  is  the  exact  theory  based  upon 

12  3 
separation  of  variables  in  elliptic  coordinates  '  '  ,  This  procedure  is  useful 

only  when  the  ratio  of  slit  width  to  wavelength  is  small  (that  is,  say,  smaller 
than  3),  Second,  there  are  approximation  procedures  that  are  valid  when  th« 
ratio  of  slit  width  to  wavelength  is  small.  Such  procedures,  based  on  the  solu- 
tion of  related  electrostatic  problems,  were  indicated  by  Rayleigh  and  have  been 

C  A  7 

developed  further  by,  e.g.,  Soramerfeld  ,  Groschwitz  and  Honl  ,  and  Bouidcamp  . 
Finally,  there  are  approximate  treatments  valuable  particularly  for  large  ratios 
of  slit  widths  to  wavelength;  among  these  are  calculations  based  on  Kirchoff 's 

procedures  and  their  various  modifications,  and  the  work  of  Schwarzschild  .  A 

9 
calculation  related  to  that  of  Schwarzschild  has  recently  been  published  by  Fox 

in  connection  with  diffraction  by  a  pulse.  It  should  be  noted  that  variational 

techniques  may  be  applied  in  conjunction  with  any  of  the  approximate  methods 

Just  mentioned. 
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The  procedure  adopted  here  is  closely  related  to  the  work  of  Schwarz- 
schild,  who  indicated  a  series  of  successive  approximations,  based  on  the 
Somraerfeld  solution  of  the  problem  of  diffraction  by  a  half -plane.  Schwarz- 
schild  considers  the  slit  configxiration  as  being  composed  of  two  diffracting 
half -planes.  Each  half -plane  responds  to  a  variety  of  excitations,  and  the 
response  in  each  case  can  be  calculated  by  use  of  the  Sonmerfeld  Green's  func- 
tion for  the  half-plane.  The  excitations  in  question  are   calculated  successively, 
as  are  the  responses.  The  first  excitation  is  the  incident  field,  and  the  cor- 
responding first  approximation  to  the  total  field  is  the  sum  of  the  fields  dif- 
fracted independently  by  each  half-plane,  according  to  Soramerfeld,  under  this 
excitation.  The  second  excitation  for  the  left  half -plane,  say,  is  the  field 
scattered  by  the  right  half -plane,  and  so  on.  This  leads  to  an  infinite  series 
of  terms,  which  Schwarzschild  shows,  is  convergent  for  every  slit  width.  How- 
ever, this  method  is  most  valuable  for  large  separation  of  the  half-planes. 
Schwarzschild  calculated  explicitly  only  the  first  term  described  above,  in  which 
the  excitation  for  each  half-plane  is  the  incident  wave  only.  This  term  will 
henceforth  be  called  the  'non-interaction  approximation'. 

The  investigation  pxirsued  in  the  present  paper  is  based  on  a  similar 
conception  of  the  action  of  a  slit,  but  the  procedure  used  here  differs  some- 
what from  that  of  Schwarzschild.  The  field  at  a  point  is  considered  as  being 
composed  of  the  incident  field  plus  a  response  field  from  each  of  the  two  half- 
planes,  or  alternatively,  as  a  geometric  optics  field  plvis  two  diffracted  fields 
(one  from  each  half-plane).  It  is  assvuued  that  each  of  these  diffracted  fields 
is  qualitatively  of  the  same  nature  as  the  known  field  diffracted  by  an  isolated 
half -plane.  In  particular,  it  is  supposed  that  the  partial  field,  i.e.,  the 
field  ascribed  to  one  of  the  half-planes,  say  the  right  half -plane,  can  be  well 
represented  in  the  region  far  enough  to  the  left  by  a  function  of  the  form 
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(e         / /r  )f(0)     (i.e.,  the  form  of  the  function  representing  the  field  dif- 
fracted by  an  isolated  half-plane)}  here  r  and  9  are  polar  coordinates  with 
pole  at  the  ed^  of  the  right  half -plane.     However,  the  angular  function  f(9), 
known  for  diffraction  by  an  isolated  half -plane,  must  be  determined  for  the  case 
of  diffraction  by  a  slit.     Now  on  the  left  half -plane,  ©  has  the  constant  value 
n.     Thus  on  the  left  half -plane,  say,  the  field  due  to  the  right  half -plane  is 
everywhere  of  the  fonr.     const.  X(e~       /  V^)i  and  this  represents,   approximately, 
the  field  of  a  line  source  situated  at  the  edge  of  the  right  half-plane,     A 
similar  remark  holds  for  the  net  excitation  of  the  right  half -plane  by  the  net 
field  ascribed  to  the  left  half-plane.     Hence  each  half -plane  must  be  thought  of 
as  being  excited  by  the  incident  field  plus  a  line  source  located  at  the  edge  of 
the  opposite  half-planej  the  strength  of  this  line  source  depends  on  the  value 
of  f(n)  for  the  opposite  half -plane.     The  total  field  can  therefore  be  expressed 
as  a  linear  combination,  with  constant  coefficients,  of  fields  of  a  well-known 
and  simple  form,  i.e.,  the  fields  which  occur  in  diffraction  of  plane-wave  and 
line-source  fields  by  half-planes.     The  coefficients  are  calculated  by  comparing, 
for  each  half -plane,  the  value  of  f(n)  as  an  excitation  and  as  a  response;  this 
yields  two  linear  algebraic  equations  for  the  two  coefficients  involved. 

As  a  resiilt  of  this  approximate  calculation,  the  diffracted  field  is 
expressed  as  the  sum  of  a   'non-interaction'  term  plus  an  'interaction'  term, 
the  non-interaction  term  being  identical  with  Schwarzschild ' s  leading  term.     Our 
interaction  term,  conveniently  expressed  in  terms  of  the  response  of  a  half-plane 
to  a  line  source  of  a  prescribed  amplitude  placed  at  the  opposite  edge,  has  a  form 
presumably  different  from  the  function  that  would  arise  when  calculating  Schwarz- 
schild 's  second  term.     This  difference  would  probably  be  a  result  of  approximating 
the  same  function  by  different  methods,   and  involves  terms  of  higher  order  in  the 
parameter,  i,e.,  the  slit  width  divided  by  wavelength.     It  seems  possible  that  the 
interaction  term  derived  here  is  more  reliable  since  it  partially  takes  into  ac- 
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count  the  total  interaction.  It  is  evident  from  the  method  used  that  the  re- 
sult derived  here  is  at  best  an  asymptotic  approximation  «md  that  it  can  be 
expected  to  be  good  when  the  ratio  of  slit  width  to  wavelength  is  large,  and 
further,  that  for  any  such  value  of  this  parameter,  it  should  be  an  improvement 
over  Schwarzschild's  first  term.  This  was  found  to  be  so  in  the  cases  for  which 
calculations  were  made* 

In  Section  II  the  relevant  known  formulas  dealing  with  diffraction  by 
a  half -plane  are  reviewed  and  developed  in  a  suitable  fona.  Section  III  is  de- 
voted to  the  calculation  described  above,  for  the  case  of  a  plane  wave  incident 
from  an  arbitrary  direction,  and  in  Section  IV  we  deal  with  the  special  case  of 
normal  incidence .  Numerical  computations  were  carried  out  for  twelve  values  of 
the  ratio  of  slit  width  to  wavelength  in  the  range  in  which  this  parameter  vairies 
between  -^1  and  '^-2,5.  This  range  substantially  repiresents  the  overlap  between 
the  upper  limit  for  which  the  exact  calculations  can  be  conveniently  handled 
and  the  lower  limit  to  which  the  asymptotic  approximation  derived  here  can  be 
pushed.  In  each  case,  the  far  transmitted  field  is  calculated  in  ten-degree 
inteirvals  for  a)  the  exact  calculation,  b)  the  non-interaction  approximation, 

and  c)  the  approximation  including  the  interaction  tenn  derived  here.  Although 

12  3 
the  calculations  for  the  exact  field  have  been  done  previously  *  *  ,  it  was 

found  desirable  to  carry  out  these  calculations  independently  using  the  Bureau 

of  Standards  Tables  relating  to  Mathieu  Functions  • 

The  n\imerical  results  are  presented  in  Tables  I  and  II,  the  field 
strengths  being  listed  in  the  former,  the  transmission  coefficients  in  the  latter. 
It  is  seen  from  Table  I  that  the  interaction  correction  derived  in  this  paper 
significantly  improves  the  non-interaction  approximation,  and  it  seems  plausible 
that  this  correction  term  may  be  used  with  confidence  for  all  slits  more  than 
one  wavelength  in  width. 

In  Section  V,  the  field  transmitted  by  a  slit  vhen  excited  by  a  line 
source  is  considered, while  Section  VT  is  devoted  to  concluding  remarks. 
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II,  T)if  fraction  by  a  single  half -plane 

In  this  section  we  review  the  well-known  solution  of  the  problem  of 

diffraction  of  a  line  source  or  plane  wave  by  a  conducting  half -plane,  and  take 

special  notice  of  certain  simple  properties  which  will  prove  useful  in  the  later 

sections.  We  consider  E  polarization,  with  an  electric  field  whose  z-component 

we  call  E  •  The  latter  must  vanish  on  the  screen, 
z 

The  field  of  a  line  source  in  the  presence  of  a  conducting  half -plane 
whose  edge  is  parallel  to  the  source  is  well-known.  Suppose  the  line  soxxrce  is 
parallel  to  the  z-axis  and  that  its  field,  in  the  absence  of  the  half -plane,  is 

(1)  E^-iH^2>(kR). 

Then  the  field  in  the  presence  of  the  half -plane  is  given  by  the  integral  ex- 
pression 


(2)  E  -  i  J  r  °  e"^  '°'^  ^  d4  -  r   e-^^l^°^^  ^  d4 


where 


r    t 


(3) 


2  /rr^ 
sinh  4^  •  — g—^  cos  ^  (©  +  9^), 

R  -  /(x-x^)  +  (y-y^)  ^ 


The  coordinates  ai^  defined  in  Fig.  1» 
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^^'oV 


i^^o'V 


P(x,y) 


Figure  1 

The  half -plane  covers  that  portion  of  the  plane  y  »  0  for  which  x  >  0, 
and  its  edge  is  coincident  with  the  z-axis.  The  source  S  is  located  at  (x^,y^) 
or,  in  polar  coordinates,  (r  ,&  ),  and  the  field  point  P(x,y)  is  denoted,  in 
polar  coordinates,  by  (r,©).  R  is  the  distance  from  the  source  S  to  the  field 
point  P|  R-  is  the  distance  from  S,  to  P,  S,  being  the  image  of  S  mirrored  in 
the  plane  of  the  half -plane. 

The  far  field  (that  is,  the  field  where  kr  »  1)  is  approximately  of 
the  fona 


(U) 


E  ^  E^^^  *  /Z    ^ 


-ikr 


2T 


yi? 


f(r^,6o,e). 


except  in  certain  regions  shortly  to  be  desciT.bed.  An  alternative  description 
is 
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Here  F(r^,e^,©)  is  e3q)ressible  in  Fresnel  integrals.  The  expression  (Ua)  is 
valid  for  all  Q  and  sufficiently  large  r,  Eqs.  (U)  and  (Ua)  can  be  derived 
from  Eq.  (2),  In  E4.  (U),  E„  ^  is  the  field  calculated  on  the  basis  of  geometric 
optics  and  is  specified  in  terms  of  three  distinct  regions,  as  illustrated  in 

Figure  2  for  a  particular  case.  Region  I  is  the  geometric  shadow  of  the  half- 

(e) 
plane,  and  in  this  region  E^^  •  0;  n  is  illuminated  by  direct  rays  only,  and 

here  E^^^  '  7T^       (kR)j  and  Region  III  is  illuminated  by  both  direct  and  re- 
flected rays,  so  that 


The  lines  separating  the  regions  are  known  as  'shadow  lines',  and  on  them 

(a) 

F    is  discontinuous.  The  angular  neighborhoods  of  these  lines  are  the  ex- 
ceptional regions  we  mentioned  in  connection  with  Equation  (U)» 


,ray  from  source 


direct  and 
reflected 

illumination 


I 


shadow  zone 


Figure  2 
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We  shall  define  E    to  be  the  diffracted  field,  everywhere  given  by 
z 


(Ub)  E   -  E^^^  ♦  eI"*^  . 


z  z 


The  second  term  on  the  right  side  of  Eq«  (U)  is  then  a  far-field  approxination 
to  E^  ,  and  since  for  kr  »  1, 

-ikr 


'^'         /S  ^  -  *  "o^''"*' 


the  far  diffracted  field  is  approximately  that  which  would  be  produced  by  an 
'inhomogeneous'  line  source  placed  at  the  edge  of  the  half-plane. 

For  the  present,  it  is  necessary  to  have  an  expression  for  the  line- 
.source  field  only  for  the  case  where  the  half -plane  is  excited  by  a  line  source 
located  on  the  extension  of  the  half-plane  (i.e.,  9  •  ti).  In  this  case  both 
the  shadow  and  reflected  zones  (Regions  I  and  III)  disappear,  so  that  E„  '^  is 
just  the  incident  field  and  f  is  given  by 

,2  Vkr  /n  sin(e/2) 


(6)      f(r^,n,6)  -  e        {~1  *  V^ 


as  can  be  seen  by  evaluating  (2)  for  large  kr.  When 


2 

exp(-inv  /2)dv 


(7)  2  y  -—  sin  I  »  1, 


/5 

'     n 

then  (6)  can  be  approximated  by 

-ikr^ 
(8)  f(r   .Ti,6)^    -     ^ 


o 
V2Tiikr     sin  » 

From  (7)  is  is  clear  that  the  approximation  (8)  is  valid  only  if  kr^  is  large 
compared  to  unity,  and  furthermore  if  the  angle  0  is  not  close  to  zero  or  2n 
(i.e.,  the  field  point  must  not  be  near  the  half-plane). 

When  the  incident  field  is  not  the  line  source  (1)  but  instead  a 
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plane  wave  of  unit  amplitude  approaching  from  the  direction  6  >i.e.,  when 

-ikr  cos(©-©  ) 
^^^  "incident  "  ®  ' 

then  the  far  field  (kr  »  1)  in  angular  regions  bounded  away  from  the  shadow 
line  is  given  approximately  by 

(10)  E  -E^^K/^  2i_  g(e^,9), 

where  the  geometrical  optics  term  E^^^     is  now  zero  in  Region  Ij  exp[^ikrcos(»-6   f] 
in  Region  IIj  and  exp[^ikrcos(9-9   )]  -  exp[2ikrcos(e+©   )J  in  Region  III.     Away 
from  the  shadow  lines,  the  angvilar  fvmction  g(9  ,6)  is  given  by 

(11)  g(o^,e)  "k  < TT ^, : 

°  '^"     I  cos  ^  (9  -  9^)       COS  J  (9  *  9^) 

Thus,  like  the  line-source  field,  the  field  produced  when  a  half -plane 
is  excited  by  a  plane  wave  can  be  expressed  as  the  sum  of  a  geometrical  optics 
term  plus  a  diffraction  term  whose  far-field  approximation  has  the  form  of  an 
inhomogeneous  line  source  located  at  the  edge  of  the  half -plane. 

III.  Diffraction  by  a  slit 

Since  an  infinite  conducting  plane  with  a  slit  is  a  pair  of  adjacent 
half -planes,  an  approximate  expression  for  the  field  transmitted  by  such  a  slit 
as  a  result  of  excitation  by  a  plane  wave  can  be  obtained  in  terms  of  the  well- 
known  results  of  the  previous  section.  As  before,  the  field  may  be  expressed 
as  the  sum  of  a  geometric  optics  term  plus  a  diffraction  term: 

(12)  E  -  E^^^  *  E^^^    . 

^  z    z     z 

In  the  three  regions  required  to  specify  El^t   shown  in  Figure  3,  the 

z 

geometrical  optics  fields  are,  respectively:   zero  inRegiorela  and  lb,  the  in- 
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cident  wave  in  Region  II,  and  the  incident  wave  plus  reflected  wave  in  Region 
Ilia  and  Illb.     The  two-half -planes  are  denoted  by  the  subscripts  a  and  b. 


Figure  3 

and  to  simplify  the  mathematics  two  coordinate  systems  are  enployed,  for  the 
time  being,  with  origins  at  the  edges  of  these  half -planes  (see  Figure  U)» 
The  plane  wave  is  incident  from  directions  ©   and  ©,   in  the  two  coordinate 
systems  respectively,  \rtiile  the  coordinates  of  the  field  point  are  (r  ,©  ) 
and  (r,  ,&,  )  respectively.  Taking  the  slit  width  to  be  2d  (d  is  what  is  common- 
ly called  the  (•  half -width  •  )>  each  edge  is  located  at  (2d,  n)   with  respect  to 
the  coordinate  system  of  the  opposite  edge* 

Tha  mathematical  formulation  of  the  problem  is  as  follows.  Given 

an  incident  field  u"^°(x,y)  »  E^  ^,  we  have  to  determine  a  field  U(x,y)  ■  E^, 

in  accordance  with  the  conditions  that 

(i)      V^  +  k^  •  0   (except  on  the  half -planes  (a)  ^d  (b))j 

(ii)      U(x,y)  -  i:u^^(x,y)  -  U^^=(x,-y)]    r"^/2exp(-ikr)f^(0),    (y  >  0)i 


-  u  - 


(iii) 

(iv) 


V{x,y)  r^  r'^^^   exp(-ikr)f„(©). 


U  ■  0     at  the  surfaces  (a)  and  (b). 


(y  <  0); 


Figure  U 

Here  r  and  ©  are  polar  coordinates  with  pole  at  the  center  of  the  slit.  For 
y  >  0  we  have  0  <  6  <  n,  while  for  y  <  0,  we  have  n  <  6  <  2jt.  Conditions  (ii) 
and  (iii)  are  to  hold  for  large  kr  uniformly  in  6  in  the  intervals  (O  <  0  <  n) 
and  (n  <  &  5  2ti)  respectively. 

Let  us  now  wirite  the  solution  in  the  form 

a    D 

The  superscript  r  denotes  'response*.  We  are  here  trying  to  express  U  in  the 

form  of  an  incident  field  plus  two  response  fields,  each  originating  in  the 

J" 

half -plane  indicated  by  its  subscript.  We  assxime  therefore  that  U„  is  regular 

in  the  exterior  of  surface  (a),  and  U.  is  regular  in  the  exterior  of  surface  (b). 


-  12  - 


r     r 
We  now  maKe  some  further  assumptions  as  to  the  nature  of  U  and  U.  ,  and  then 

try  to  determine  these  functions  approximately.  The  success  of  this  procedure 

will  then  depend  on  the  degree  to  which  the  regularity  and  boundary  conditions 

are  fulfilled  by  the  resulting  approximation  to  the  function  U. 

r   r 

We  assume  that  the  response  fields  U  ,  U,  have  respectively  the  same 

behavior  as  the  known  response  which  half-planes  (a)  and  (b)  would  make  to  the 
incident  field  if  each  were  isolated.  To  express  this  assumption  mathematically 
we  first  define  functions  E  ^  and  E   by  the  following  equations  (cf.  Fig.  3): 

U^  .  -  -U^°(x,-y)  +  Ef  >^^  in  Ilia,  Illbj 

o^  D  Z 

(12a)         U^  ,  -  -U^"^(x,y)  *  E^^*^  in  la,  lb; 

S  •  D  Z 

U^   .  E'^a^db  otherwise. 

a,b    z 

r   da 
The  above  equations  are  to  be  i^ad  twice,  once  as  applying  to  U  ,  E   and 

8.    Z 

regions  Ula  and  la,  and  again  as  applying  to  uf ,  E„  and  regions  Illb  and  lb, 

D    Z 

Then  the  assumption  of  a  response  behavior  analogous  to  the  isolated  half -plane 
responses  is  seen  to  be  equivalent  to  the  postulation  of  the  asymptotic  formulas 

(12b)         Sf '"^^  (^a,b)"'^'  «^(-i^a,b)^a,b(»a>^b)- 

These  formulas  are  assumed  to  be  vaHd  for  large  kr  or  kr^  respectively,  in 

da 
angular  regions  bounded  away  from  the  shadow  lines.  Thus  for  E^  ,  Equation  (12b) 

is  assumed  to  be  valid  for  9     bounded  away  from  the  lines  labeled  a  in  Figure  3. 

3. 

For  E*  the  lines  labeled  B  are  excluded.  The  behavior  implied  in  Equations 
z  '^ 

(12a)  and  (12b)  actually  occurs  in  the  case  of  a  single  half -plane.  Of  course 
we  are  not  assuming  that  the  functions  F^(©^)  and  F^(&^j)  are  the  same  as  those 
for  an  isolated  half -plane. 

The  assumptions  we  have  made  thus  far  can  be  siimmarized  by  saying  that 
we  have  assumed  that  the  function  U  defined  in  our  mathematical  formulation  can 
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be  represented  in  the  form  U»U         +U     +Ur,  where  the  fiinctions  U     and  U, 

a    D'  a     D 

are  restricted  by  the  requirements  expressed  in  (12a)  and  (12b)  and  by  the 
requirements  that  U  ,  be  regular  everywhere  except  on  the  surfaces  (a)  and 
(b)  respectively. 

It  should  be  remarked  here  that  the  crux  of  our  argument,  and  the 
approximation  involved  in  it,  has  not  yet  appeared.  In  fact,  the  physically 
plausible  assumptions  mad©  so  far  can  actually  be  proven  rigorously}  they  are 
mathematical  consequences  of  the  conditions  of  the  problem.  Although  it  is 
sufficient  to  make  these  physical  assumptions  outright,  since  when  we  obtain 

our  approximate  solution  we  shall  verify  them  in  any  case,  it  seems  useful  to 

inc 
indicate  a  proof  here.  This  would  consist  of  expressing  the  response  to  U   , 

by  means  of  Green's  theorem,  in  the  form  of  two  integrals  whose  surfaces  of 
integration  are  (a)  and  (b)  respectively.  Each  of  these  integrals  is  a  regular 
wave  function  outside  its  surface  of  integration.  The  integrals  represent  dis- 
tributions of  line  sources.  Since  the  density  of  distribution,  or  current 
density,  can  be  expressed  in  terms  of  the  normal  derivative  of  U,  its  asymptotic 
behavior  can  be  shown  [using  the  conditions  (i)  -  (iv)]  to  approach  a  value  equal 
to  twice  the  normal  derivative  of  the  incident  field.  From  this  one  can  show 
that  the  surface  integrals  we  have  described  behave  for  large  kr  ,  kr,  in  the 
way  we  have  postulated  for  our  functions  U^  and  uf.  Thus  the  Green's  function 
integrals  actually  provide  a  representation  of  the  exact  solution  in  the  form 
we  have  assumed. 

We  now  proceed  with  the  essentially  approximate  part  of  our  argument. 
From  equation  (12a)  we  know  that  near  half -plane  (a),  the  function  U^^  is  given 

by  E*^^.  The  latter  function,  in  turn,  is  now  seen  to  be  equivalent  (near  half- 
z 

(2) 
plane  (a))  to  some  constant  multiple  of  a  line  source  H^  '(kr,  ),  provided  the 

slit  is  large  compared  to  wavelength.  This  follows  from  equation  (12b)  if  we 

note  that  the  function  r"    exp(-ikr,  )  is  a  constant  multiple  of  the  function 
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H   (kr,  ),  for  large  kr,  ,  and  also  that  along  half-plane  (a)  we  have  6-.  ■  « 

and  F.  (6.  )  ■  F,  (jt)  ■  constant.  We  now  observe  that  the  meaning  of  the  above 
b  b     D 

result  is  that  hsLLf-plane  (a)  is  being  excited,  in  effect,  by  the  incident 
field  and,  in  addition,  by  a  line  source  of  undetermined  strength  located  at 
the  edge  of  half -plane  (b).  We  argue  that  the  response  of  (a)  must  be  equal 
to  the  known  response  of  an  isolated  half-plane  to  the  incident  field  plus 
some  as  yet  undetermined  multiple  of  the  known  response  which  such  a  half- 
plane  makes  to  a  unit  line  source.  Similarly  we  find  that  half -plane  (b) 
is  being  excited  by  a  line  source  located  at  the  edge  of  (a),  in  addition 
to  the  given  incident  excitation.  We  are  therefore  led  to  the  conclusion 


.r 


res- 


that  U,  can  similarly  be  expressed  as  a  linear  coinbination  of  the 
ponse  to  the  incident  field  and  the  response  to  a  line  source.  Let  U^  be 
the  known  response  of  half-plane  (a)  to  the  given  incident  field  when  half- 
plane  (b)  is  absent,  and  let  IT.  be  the  known  response  of  the  half-plane  (a) 
to  a  line  source  w^  H^  (kr,  )  again  when  (b)  is  absent.  This  line  source  is 

located  at  r  ■2d,  ©  ■  n  in  the  r  ,  0  coordinate  system.  We  have  seen  that 
a     '  a  a'  a 

uf  has  the  form  const,  X^  H^  ^^b^  ^^   points  on  half -plane  (a).  Let  C^  be 
the  unknown  constant  involved.  Also  introduce  the  known  functions  Uv^»  ^ha* 
and  the  unknown  constant  C.  in  a  similar  manner:  U,  .  is  the  response  of  iso- 
lated half-plane  (b)  to  the  incident  field^uf  is  the  response  to  the  excita- 
tion ir  H^  (kr  ),  and  C.  is  the  constant  involved  in  the  asymptotic  repre- 
2i  o     a       D 

sentation  of  U^  at  points  on  half -plane  (b).  Then  %re  may  write 
a 

r  j^   r       r 

"I  ■  "ai  ^  ^a  ^ab  ' 

(")         "b  ■  <i  *  Vba  . 


U  -^  K 


inc    ai    bi 


>    *  ^a  <b  *  ^b  "ba  • 
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<'^^'v^ 

g(»ao»®a^  *     V^2^ 

f(^bo>®b)  *     S^(^^ 

This  representation  holds  at  any  point  in  space.*  To  complete  the  specification 

of  U,  U  and  U,  we  must  detenrdne  the  unknown  constants  C  and  C.  ,  Before 
'  a     b  a     D 

turning  to  this  problem  we  consider  the  quantities  E  ,  E   given  in  Eqs.  (12b). 
These  equations  may  now  be  rewritten  in  greater  detail  in  view  of  Eqs.  (13) »  if 
WB  use  the  asymptotic  forms  of  the  diffracted  fields  of  the  half -plane  for  plane- 
wave  and  line-source  excitation.  Using  the  notations  of  Section  II,  Eqs.  (10) 
and  (11)«  and  Eqs.  (U),  (6),  (7)  and  (8),  we  find 

-ikr 
(13a)         Ef^  ^  )/%    1^  rg(».o>6.)  *     C,f(2d,n,9j]  , 

(13b)         Ef  ^  /^  ^~;I  |g(^bo>®b)  *  C^f(2d,n,e^)l  . 

These  formulas  hold  if  we  let  r  or  r,  tend  to  infinity  while  ©^  (or  ©^)  are 
bounded  away  from  the  lines  a  (or  the  lines  p)  shown  in  Figure  3. 

Let  us  also  note  at  this  point  a  significant  regrouping  of  Equations 
(13),  (13a),  and  (13b).  We  can  write 

(13c)  u  5  E,  •  E,(non)  +  E,(int) 

with 

(13d)         E^(non)  -  "i^c  *  ^  "  ^i 

and 

(13e)         E^(int)  -  cXb  *  %^la  ' 

The  non-interaction  terra  E     '(non)  is  exactly  the  first  term  of  the  Schwarzschild 

z 

series  mentioned  in  Section  I,  and  represents  the  excitation  of  each  half-plane 
by  the  original  plane  wave  only.     It  may  be   a  good  approximation  for  extremely 


«     It  is  to  be  noted  that  it  is  an  approximate  representation  in  terms  of 
known  functions. 
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wide   slits,  in  which  case  the  half -planes  have  little  effect  on  each  other.     For 
narrower  slits,  where  the  interaction  of  the  half -planes  cannot  be  neglected, 
the  correction  terms  El     (int)   is  added. 

If  we  take  oixr  field  point  to  infinity  in  the  lower  haJJ'-space  in  a 

fixed  direction  not  parallel  to  the  shadow  lines  a  and  p,  then  we  enter  either 

Y»  inc         dfl  r         db 

Region  I     or  L  .     In  Region  I    we  have  IT   •  -U         +  E       and  U,    ■  E     ,  while  in 

Region  L    we  have  uj"  »  -U^"°  ♦  E^^  and  U^  -  E^^   [cf.  Eq.   (12a)].     Thus  in  either 

D  D  Z         £1     Z 

case  we  find  that  U  ■  E  is  given  by  the  sum  of  E_  and  E„  ,  which  are  a  pair 

Z  Z        z 

of  functions  of  the  form  indicated  in  (12b) j  these  functions  are  now  to  be  further 
specified.  Furthermore  in  Regions  I  and  L  the  geometrical  optics  field  as  we 

el        D 

specified  it  in  connection  with  Fig.  3  is  zero,  and  hence  the  total  field  is  also 
the  diffracted  field  as  defined  in  Eq.  (12),  Consequently,  if  we  pass  to  infinity 
in  the  lower  half -space,  in  a  direction  not  parallel  to  a  shadow  line,  the  follow- 
ing formulas  hold  in  virtue  of  Equations (13a)  -  (13e): 

(lU)         E  s  E^**^  -  E!j'^^(non)  ♦  E^^^^int), 


where 

(15a)         E^^^non)  -  /g  . 


and 


r  -ikr  -ikr^j 


-ikr  -ikr, 

a  .    b 


(15b)         E^'^^int)  a  /^  JS-----  C/(2d,n,©^)+  1—-  C^^f  (2d,n,©^,)  I  . 

We  now  turn  to  the  determination  of  the  constants  C  and  C,  •  We  have 
been  regarding  half-plane  (a)  as  subject  to  the  excitation  of  the  incident  field, 
and,  in  addition,  to  excitation  by  a  fictitious  line  source  located  at  the  edge 
of  (b).  This  line  source  gives  the  effect  of  the  function  U^.  Equation  (13b) 
shows  that  near  half -plane  (a)  we  have  U^  s  E^  j  an  expression  for  the  latter 
is  given  in  (I3b).  Since  (n/2ikr^)-'-/^exp(-ikr)  '^^  (n/2i)H^^^kr^),  we  see  that 


-  17  - 


e'^ 'o  constant  X  (n/2i)H^^^(kr.  ). 

2  O      D 

Ths  constant  of  propoirtionality  is  (from  (13b)) 

(16)  g(»bo»"^  *  C^f(2d,n,ji). 

Hence  the  partial  far  diffracted  field  produced  by  the  half -plane  (a)  in 

specific  response  to  the  excitation  by  the  fictitious  line  source  at  the 

edge  of  (b)  is  equal  to  the  constant  (16)  times  the  well-known  response 

[cf ,  (U)]  to  a  homogeneous  line  source,  and  is  given  by 

-ikr 

(17)  [^(©^^,it)  +  C^^f (2d,n,n)l /^ f(2d,Tt,9^). 

'   a 

Now  from  Fq.  (13)  we  find  a  different  expression  for  the  partial  field  (17), 
namely  ^^^ 

(18)  C„/X  ^—J-    f(2d,n,a  ). 


a 


^ 


'   a 


Equating  expressions  (17)  and  (18),  we  have 

(19)  C^  -  g(©bo,n)  +  C^jf(2d,n,n). 

Proceeding  as  abo^i•  with  respect  to  the  diffracted  field  produced  by  half-plane 
(b)  due  to  excitation  by  the  fictitious  line  source  at  the  edge  of  (a),  a  similar 
relation  is  obtained: 

(19a)         C^  -  g(»ao»''^  *  C/(2d,n,n). 

Equations   (19)     and  (l9a)  are  two  linear  inhoraogeneous  equations  in 

the  desired  coefficients  and  yield 

g(»^o,")  +  g(6^Q,ii)f(2d,n,n) 

C  •  A 

^         1  -  [f(2d,n,r)]^ 

^^°^  g(»ao»"^  *  g(©^^,n)f(2d,n,n) 

C  •  -5 . 

^         1  -  [f(2d,n,n)]'^ 
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The  final  formula  for  the  field  at  ary  point  is  given,  in  terms  of 
line-source  and  plane-wave  fields,  by  equations  (13)   and  (20 )• 

In  (20)  we  may  substitute  the  approximation  (8)   for  f(2d,n,n)  if  the 
slit  width  is  not  less  than  a  wavelength,  and  then,  substituting  the  expression 
(11)  for  g(&  ,e)  we  have  finally 

/  ^  ^-2ikd  \    /  ^-iUkd\-l 

^a-    ^-  -in^a^^/.;  *  ,3in(»^/2)  A^  J  ^  '  "  ^^^^ 
(21)  ^ 


C--     l-nsin(9^y2;\^.^(^^^/2)VCTi)     (^"^^^^^ 

The  field  representation  which  includes  the  interaction  tern  should 
be  a  better  approximation  than  the  'non-interaction •  representation,  and  hence 
may  be  expected  to  hold  for  smaller  slit  widths  than  the  latter.  It  also  seems 
to  have  an  advantage  over  the  well-known  solution  in  the  form  of  an  infinite 
series  of  Mathieu  functions,  in  that  it  is  an  intuitively  clear  representation 
in  terms  of  known  elementary  functions. 

For  convenience  of  expression  the  results  derived  in  this  section  have 
been  given  in  terms  of  a  double  set  of  coordinates.  For  computation,  however, 
it  is  more  convenient  to  express  these  results  in  a  single  coordinate  system 
with  origin  midway  between  the  half -planes  (see  Fig.  5). 

Because  the  equations  of  transformation  between  the  various  systems 
of  angular  coordinates  described  above  are  different  for  y  >  0  and  y  <  0,  the 
case  where  the  field  point  and  the  incident  wave  are  on  the  same  side  of  the 
plane  y  -  0  cannot  be  considered  together  with  the  case  where  the  field  point 
and  incident  wave  are  on  opposite  sides.  Only  the  latter  case  will  be  consider- 
ed here,  and  will  yield  the  field  transmitted  by  the  slit. 
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+d 


Figure  $.     Single  coordinate 
system .     All  angles  are  the 
acute  angles  measured  from 
the  nearest  normcil  and  are 
considered  positive   if  measured 
coianter-clockvriLse  with  respect 
to  this  normal,   and  negative 
if  clockwise. 


Let  the  incoming  plane  wave  be  incident  from  above  the  slit  and  the 
field  point  P  be  below*     Further,  let  0     represent  the  acute  angle  between  the 
normal  to  the  incoming  plane  wave  and  the  normal  to  the  plane  of  the   screen 
(measured  from  the  positive  y-axis),  and  let   (r,5>)  represent  the  field  point  P, 
where  ^  is  the  acute  angle  between  the  field  point  and  the  nonnal  to  the  screen 
(measured  from  the  negative  y-axis) »     All  angles  are  considered  positive  if 
measured  counterclockwise  with  respect  to  the  normal  and  negative  if  clockwise. 
The  relations  between  these  coordinates  and  those  in  each  of  the  systems  pre- 
viously used  is,  in  general,   quite  complicated,  but  when  P  is  far  from  the 
edges  as  compared  to  the  slit  width,  i.e.,  when 

approximate  relations  between  them  can  be  simply  stated.     Remembering  the  con- 
vention that  angles  measured  counter-clockwise  are  positive  and  angles  measured 
clockwise  are  negative,  we  have 


a)     r    ai  r  +  d  sin  ^ 

(22)  c)     © ^  1  n  -  (21 

a       c 


b)       r,  i«i  r  -  d  sin  (if 
d)       e^  ^  I  n  +  (2 


f)       9 


bo  ■  ? 


*     ^r 


-  20  - 


The  difference  between  (c),   (d)   and  (e),   (f)  is  due  to  the  fact  that 
(f  is  measured  from  the  negative  y-axis,  vrtiile  (f     is  measured  from  the  positive 


y-axis» 


Substituting  the  relations  (22)  into  (11)  yields 


aC   a 


K»».»&J-  -  h\  — r *  — r ^ ' 

^"    t^in  J  (i?-^^)  cos  ^  {9^<?^)  ^ 


(23) 


g(®bo'®b^"-^< 


7^ 


-1 

T 


+    T 

^sin  J  (0-^^)         cos  J  (P+0q) 


Similarly,   substituting  (22)  into  (6)  yields 

p2  vW"    [cos(gi/2)+sin(9'/2)] 


f(2d,n,©  )-e*^2*^'^«^"^  J 


-1+  y/2r 


V. 


J 


exp 


(2U) 


f(2d,n,9^)*e-"^^^^<'^-l*/2r 


■2  y^dTn    [cos(9r/2)-sin(9l/2)] 


exp 


0 


[^inv2/2J 
[.inv2/2] 


dv 


For  the  values  of  0  where  the  appro xLraation  (8)  is  valid,   i.e.,   away  from  the 
shadow  lines  6     "0  and  ©.    ■  0,  these  expressions  reduce  to 


-i2kd 


(25) 


f(2d,Tl,»^)  ^  - 


f(2d,n,©^)  ^  - 


e      1 

/^nikd  cos  ^  +  sin  i- 


-i2kd 


v^2aikd  cos  2-  -  sin  ^ 


Finally,  the  coefficients  C^  ^^  ^^  become 


-  21  - 


( 

/?  /n 


T       e 

"  Unikd 


^-i2kd 


cos (0/2)+  sinCgr  /2)       cos(0  /2)-  sin(0  /2)     y^l^Jira 


^^^^  ^  f  ^ 

P  ^/"     .  .  1  1  e-i2^d 

1  -  ^"_     1     cos(;2^/2)-  sin(0  /2)       cos(0  /2)+  sin(9>  /2)    v^I^IiEd 

Actually  we  have  not  verified  that  oxir  field  has  the  required 
asymptotic  behavior  for  large  kr    (as  specified  in  the  conditions  (ii)   and  (iii)^, 
i.e.,  that  in  terms  of  the  coordinates  of  Figure   5  we  have,  for  r  -c>oo, 

E     -E^"*'  -  E^*^®<'*'®'^-(nkr)-^/2exp(-ikr)F(i2),  y  >  0   , 

z        z 

and 

E     ^  (nkr)"-'-/^  exp(-ikr)F(9f),  y  <  0 

z 

uniformly  in  Qf,  To  verify  this  behavior  we  recall  the  representation  given 
in  (13c)  -  (13e).  In  those  equations  E_(non)  is  the  leading  term  of  Schwarz- 
schild's  solution.  The  asymptotic  form  of  this  term  causes  some  difficulty  at 
g(  .  "  for  a  reason  shortly  to  be  indicated.  However, Schwarzschild  showed  that  it 
has  the  behavior  we  require  [he  did  this  in  the  course  of  obtaining  an  equation 
equivalent  to  our  equation  (27)3 •  ^^  interested  reader  is  therefore  referred 

to  ref.  8  for  a  proof.  Now  the  cause  of  the  difficulty  mentioned  above  is  that 
for  plane-wave  excitation  the  cylindrical  wave  representation  [cf ,  Eq.  (10)] 
for  the  diffracted  fields  breaks  down  near  shadow  lines,  and  the  function 
g(©  ,©)  appearing  in  Eqs,  (11),  (13a)  and  (13b)  tends  to  infinity  [cf,  the 
restriction  as  to  direction  given  after  Eqs,  (13a)  and  (13b)].  On  the  other 
hand,  no  shadow-line  phenomenon  occurs  in  connection  with  the  line-source 
response,  so  that  the  terra  E  (int)  in  Eqs.  (13a)  and  (I3b)  can  be  readily 
shown  to  have  the  proper  asymptotic  behavior,  by  simply  employing  the  rela- 
tions (22).  In  the  next  section  the  relevant  calculation  is  carried  out  for 
the  case  of  a  normally  incident  wave. 
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It  seems  worthwhile  to  remark  that  in  treating  E  (non),  if  we  pass 

z 

to  the  limit  r  -i>oo  with  (2  f   ^  ,  thus  avoiding  the  shadow-line  direction,  and 
then  let  J?  ->  ^  ,  the  latter  limit  is  readily  found  to  exist.  This  is  a  con- 
sequence of  the  fact  that  under  these  circumstances  &  and  ©,  approach  the 
shadow  line  from  lit  and  unlit  sides  respectively,  so  that  the  divergent  parts 
are  cancelled.  The  situation  is  illustrated  in  Eq.  (27)  of  the  next  section. 


IV.  Numerical  results  on  the  normally  incident  plane  wave 

To  apply  the  results  of  the  previous  section  to  the  special  case  of 
a  plane  wave  normally  incident,  it  is  necessary  only  to  set  0^  =  0.  Making 
this  substitutiong  we  find  from  (23),  (22a),  (22b)  and  (l5a)  that  the  non- 
interaction term  of  the  far  transmitted  field  for  this  case  is,  for  r  ->oo, 
with  gr  y  J  , 


E^  '(non)  -  V  ^ 


-ikr    -ikd  sin^J 


e 

TIT 


yE*   I  y  sin  w-   cos 


(27) 


+ikd  sin^ 
' "^ 


which  can  be  written  in  the  somewhat  simpler  form 

.   .  -ikr  sin(kd  sin^  -  &)  +  i  sin(kd  sin^  +  V) 

(2«)      E^^^non)  «  1_  -l^ . 

ynkr 

Further,  if  we  set  0  -  0  in  (26)  the  coefficients  simplify  to 

-  y/?  /n 

1  +  * 


/Unikc! 

Now  we  substitute  (29)  and  (2U)  into  (I5b)  and  obtain  the  interaction  contribu- 
tion to  the  far  transmitted  field: 
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E(^'(lnt).-l-_.l_j5^< 


^+ikdsinp 


v/nkr 


/Uidkd 


-1+  y^ 


(30) 


-ikdsin0 


+  e 


-1+  v^ 


.2  v/H/lT  [cos(^/2)+sin(9(/2)] 
exp(-in  V  /2)dv 


.2  VOT^   [cos(9(/2)-Ein(p/2)] 

2 
exp(-iji  V  /2)dv 


When  (25)  is  applicable,  i.e.,  avisy  from  the  lines  (2  »  ±  i,  this  equation  reduces 
to 


(31)  E^^\int)       ® 


-ikr 


V'nkr 


/iinikd  +  e 


2  -i2kd         cos(kdsin^  -  W-icos(kdsirV  + 


cos  0 


r  • 


The  above  equations  enable  us  to  calculate  the  function  F(0)  defined  by  the  equation 
(31a)     E^  «  E**  ->'  (nkr)"-'-/^  exp(-ilcr)F(5(),         (r  ^oo,  y  <  O). 

Let  us  now  calculate  F(^)  in  the  limit  {?  ->  0.  Using  (28)  and  (31)  ve  find 


(31b)     F(0)  S  i  <  2kd  -  1  +  i(2kd  +  1) 


'  + 


2(l-i)exp(-2ikd) 


/Unikd  +  exp(-2ikd) 

Use  of  equation  (30)  yields  a  somewhat  more  accurate  result.  For  wave  functions 
vanishing  on  the  screen  it  is  well  known  that  there  is  a  simple  relation  between 


F(0)  and  the  transmission  coefficient  T.  One  can  show  that 
(31c) 


r"/2    2 


,  _  Re(l-i)F(0)  .  'io     1^1  '^ 


Since  our  function  u  is  only  approximately  correct  and  only  approximately  satisfies 
the  boundary  conditions,  we  have  calculated  T  by  both  methods  indicated  in  equation 
(31c). 

To  determine  the  accuracy  of  the  approximation  procedure  developed  in 
tSiis  paper  we  shall  use  equations  (lU),  (28),  (29)  and  (31)  to  derive  numerical 
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values  for  the  field  strengths  and  then  compare  the  resiilts  with  the  values  com- 
puted from  the  exact  solutions  in  terms  of  ^lathieu  function.  From  the  exact 

solution,  the  far  diffracted  field  transmitted  by  a  slit  when  excited  by  a  plane 

12  3 
wave  normally  incident  can  be  written  as  a  series  in  Mathieu  functions  '   *  : 

(32)     li^  E^  -  ^   V^  E  r^'^f   A2^.r^2r.l(^'^''^  '  7^\  » 
T~><x>  yvkr     I  I 

>diere  «         /  o  t  \ 

(33)         Aj^^i- — :; — — ^ • 

^  '   ^^o,2r+l 

The  notation  used  here  is  that  employed  in  the  Bureau  of  Standards  Tables^  and 

12 
the  definitions  agree  with  those  of  Stratton  .  Equations  (32)  and  (33)  have 

3 
been  derived  in  a  straightforward  manner  fol] owing  the  procedures  of  Skavlem  , 

except  for  a  difference  in  notation.  The  far-field  form  (32)  is  obtained  by 

letting  the  elliptic  coordinate  corresponding  to  r  go  to  infinity  and  factoring 

out  the  common  asymptotic  forni  of  the  functions  of  that  coordinate. 

The  usefulness  of  the  exact  solution  (32),  (33)  is  limited  to  slit  widths 

where  kd  <  10,  since  from  very  large  slits  too  many  terms  in  the  expansion  (32) 

must  be  considered.  On  the  other  hand,  the  approximation  (28),  (30)  is  valid 

only  for  wide  slits  (kd  large),  where  each  edge  is  in  the  far  field  of  the 

opposite  edge.  Fortunately,  there  exists  a  region  of  overlap  (3  <  kd  <  10)  in 

which  the  approximate  solution  can  be  compared  with  the  exact  one.  In  this 

range  of  slit  sizes,  twelve  values  were  chosen  for  which  the  coefficients  of 

the  Mathieu  functions  are  tabulated  in  the  Bureau  of  Standards  tables,  so  that 


*  The  Mathieu  functions  themselves  are  not  tabulated  as  functions  of  the  arguments, 

( pT**n  ^  p  ? 
but  only  the  coefficients  DOpk  i  (^  <^  )  ^^°"^  which  the  Mathieu  functions  are  de- 
termined by  the  formula 

So^  -L(kV,x)  =  XI  Do(|^*i^k^d^)sin(2k+l)x. 

kaO 
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interpolation  was  avoided.  The  various  quantities  appearing  in  the  fonnilas 
(32)  and  (33)  were  taken  correct  to  five  significant  figures  and  the  intermediate 
results  rounded  out  to  fourj  the  far-field  amplitude  F((Z)  is  given  to  three  de- 
cimal places,  A  spot  check  reveals,  in  some  instances,  errors  of  one  or  two 
units  in  the  last  decimal  place}  this  is  due  to  the  large  nunber  of  operations 
involved , 

In  order  to  check  the  accuracy  of  our  approximate  solution,  the  far- 
field  amplitude  F(9')  was  calculated  in  accordance  with  our  approximation.  For 
an  understanding  of  the  tables  subsequently  to  be  presented  one  must  bear  in 
mind  that  these  approximate  computations  are  intended  only  for  purposes  of  com- 
parison with  the  exact  solution,  and  that  discrepancies  are  in  most  cases  apparent 
already  in  the  second  decimal  place.  It  should  also  be  noted  that  our  formula 
(30)  involves  the  Fresnel  integral.  For  large  enough  kd  this  integral  can  be 
approximated  in  terms  of  elementary  functions,  provided  the  angle  9'  is  not  too 
close  to  90  .  In  order  to  control  this  behavior,  we  started  our  computations 
with  the  angle  5^  =  0,  and,  using  the  elementary  approximation  to  the  Fresnel 
integral,  we  proceeded  by  10  increments,  comparing  each  result  with  the  Mathieu- 
f unction  result.  When  our  level  of  accuracy  decreased  sharply  we  reverted  to  the 
use  of  the  exact  Fresnel  integral,  working,  however,  only  to  two  decimal  places. 
This  was  found  adequate  to  reduce  the  discrepancy  to  the  same  level  as  that  ob- 
tained for  smaller  angular  values  of  0,  For  ^  '  90     the  result  to  two  decimal 
places  was  F{^)   =  0  (except  for  the  smsillest  kd  we  used);  this  is  in  agreement 
with  the  boundary  conditions. 

Table  I  lists  the  values  of  the  far  field  at  10  intervals  from  {?  =  0 
to  (j(  -   90°  for:  a)  the  exact  solution  (32),  (33);  b)  the  non-interaction  approx- 
imation (28)}  and  c)  the  approximation  corrected  for  the  interaction  teriu  (30). 
All  field  strengths  are  divided  by  the  radial  factor  e    /  /nkr  ,  It  is  seen 
from  Table  I  that  the  interaction  terni  can  by  no  means  be  neglected  in  this 
range  of  slit  sizes. 
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In  addition  to  the  calciilation  of  the  far  Tields  we  also  computed  the 
transmission  coefficient 


(31i) 


T  -    I  \F\^dijl    /kd. 


For  the  exact  F{^)   this  quantity  is  also  given  by  the  formula 

(3Ua)  T  -  Re(l-i)F(0)/2kd, 

but  the  above  formulas  need  not  necessarily  agree  if  F{(jS)   is  incorrect.  In  ovir 
computations  we  calculated  approximate  values  of  T  by  the  trapezoidal  rule  from 
the  approximate  F(!Z)  of  Table  I,  using  equation  (3ii )  .  To  control  the  intrinsic 
error  due  to  numerical  integration  we  also  performed  the  same  calculation  using 
the  exact  values  of  F(^)  from  Table  I,  and  compared  these  with  the  values  obtained 
from  Table  I  via  equation  (3ii)  .  The  discrepancy  was  not  more  than  ,001.  (Use 
of  Simpson's  rule  instead  of  the  trapezoidal  rule  was  found  to  produce  greater 
errors  in  some  cases.)  In  addition  to  these  results,  it  seemed  of  interest  to 
compile  the  values  which  are  obtainable  using  our  approximate  solution  in  con- 
junction with  formula  (3Ua),  In  this  connection  F(0)  may  be  obtained  using 
equation  (31).  This  leads  to  equation  (31b),  and  so  to  the  formula 

(3Ub)  T  -  1  ^    Rer2(l-i)^exp(-2ikd)1[l  +  ^^P^^^^H      [2kdVirneJ      . 

The  approximate  value  of  F(0)* involved  here  is  obtained  by  using  an  asymptotic 
evaluation  of  the  Fresnel  integral  occurring  in  (30),  and  it  should  therefore  be 
compared  with  the  values  obtained  by  numerical  integration  from  Table  I.  A  more 
accurate  reflection  of  our  approximate  solution  is  obtained  by  using  equation  (30) 
itself  in  conjunction  with  (3ija).  This  gives  .he  approximate  formula 
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00 


(3a.c)  T  - 1  *  ^ 


Re  Uexp(-in/U)  exp( ^)dv 

1       "  1-2  ^/WR  J 


1  +   (Ujikd)'-^/^  expPi(2kd+  J)l 


vrtiich  turns  out  to  be  in  error  by  (-0.005)   for  kd  =  3,   i.e.,  when  the  slit  width 
2d  is  3/n  times  the  wavelength  X.     In  contrast,  the  simpler  formula  (3ijb)  is  in 
error  by  (-0.012)  for  this  value  of  kd. 

Of  special  interest  is  the  asymptotic  formula  obtainable  when  we  retain 
only  the  leading  term  of  (3Ub),     We  find 

(5Ud)  T(asymptotic)   =  1  =  ^-r^  cos(2kd-n/U). 

/JT  (kd)^/2 

This  simple  result  turns  out  to  be  in  error  by  (-0.02U)  at  kd  =•  3.  For  kd  =  8 
however,  the  error  is  (-0,001)  which  should  be  compared  with  the  error  of  (+0,001) 
resulting  from  all  the  previous  formulas.  All  the  qu;nntities  we  have  described 
here  sre  given  in  Table  II,  where  we  have  also  included  the  values  given  by 
Skavlem  for  3  <  kd  <  8. 

It  should  be  noted  that  the  accuracy  of  the  approximate  transmission 
coefficient  exceeds  the  accuracy  of  the  far-field  pattern.  This  is  due  to  the 
fact  that  in  computing  the  transmission  coefficient  from  the  forward  amplitude 
the  error  in  the  latter  is  divided  by  2kd,  A  better  conception  of  the  accuracy 
is  fonned  if  one  considers  the  percentage  eiror  of  the  deviation  from  xmity. 
For  kd  =  3  the  percentage  error  of  the  deviation  from  unity  of  (3Uc)>  (3^b)  and 
(3Ud)  are,  respectively,  17%,  Ul%  and  82%,  For  kd  -  8  all  the  formulas  differ 
by  about  ,001  from  the  true  value  of  1,023,  Thus  the  percentage  error  is  of  the 
order  1/23  or  U%.  Actually,  if  we  refrain  from  rounding  off,  we  obtain, instead, 
percentage  errors  of  1%,  2 /o  and  U/o  for  kd  =  8   (cf.  Table  II). 

It  seems  useful  to  present  at  this  point  a  brief  analysis  of  the  sig- 
nificance of  the  successive  simplifications  which  lead  from  equation  (3'4c), 
through  (3itb),  to  equation  (3Ud),  It  will  be  recalled  that  we  take  account 
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of  interaction  in  our  theory  by  regarding  each  half -plane  as  responding  to  a 
secondary  excitation.  This  secondary  excitation  is  taken  to  be  a  line  source, 
and  the  strength  of  the  line  source  is  then  itself  determined  by  taking  account 
of  the  interaction.  The  secondary  response  of  each  half-plane  is  thus  a  product 
of  two  factors;  it  is  the  response  to  a  unit  line  source  multiplied  by  an  assigned 
source  strength.  One  may  now  approximate  each  of  these  factors.  First  of  all  we 
can  ignore  interaction  in  determining  the  strength  of  the  line  source*  The  far- 
field  form  of  the  response  of  each  half-plane  to  the  incident  wave  itself  has 
the  form  of  a  line-source  field.  This  line-source  field  is  then  regarded  as  the 
secondary  excitation  for  tte  other  half -plane.  An  additional  simplification  is 
possible  in  calculating  the  far-field  form  of  the  secondary  response  to  a  line 
source.  It  will  be  shown  that  this  is  effecti\'ely  the  same  as  the  response  to 

a  plane  wave.  In  fact,  the  far-field  angular  amplitude  in  a  certain  direction  0, 
say,  can  be  expressed,  by  reciprocity,  as  the  response  to  a  plane  wave  coming 
from  the  direction  ©.  The  latter  response  is  to  be  evaluated  at  the  source  point. 
This  statement  is  strictly  correct  so  far.  But  when  the  source  point  is  far  firom 
the  edge  of  the  excited  half-plane,  we  now  have  available  the  simplification  in- 
volved in  replacing  the  plane-wave  response  at  the  source  point  by  its  far-field 
form  (provided  the  direction  9  is  not  too  close  to  the  responding  half-plane,  in 
which  case  the  horizontal  direction  would  be  near  a  shadow  line).  We  will  then 
be  calculating  the  far-field  response  in  the  horizontal  direction  due  to  an  in- 
cident plane  wave.  This  can  be  expressed,  again  by  reciprocity,  in  terms  of  the 
far-field  response  of  the  half -plane,  in  the  direction  ©,  due  to  excitation  by  a 
plane  wave  traveling  horizontally,  i,e.,  coming  from  the  direction  of  the 
opposite  edge.  Thus  in  the  far  field  the  response  of  a  half -plane  to  a  distant 
line  source  is  effectively  the  response  to  a  plane  wavs  of  suitable  strength, 
despite  the  fact  that  the  excitations  differ  markedly  along  the  half -plane  except 
in  the  neighborhood  of  the  edge,  and  the  responses  differ  markedly  except  in  the 
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far-field.  This  phenomenon  is  presumably  connected  with  the  fact  that  for  large 
kr  and  kd  the  wave  number  k  is  large  in  effect,  so  that  geometric-optical  con- 
siderations are  valid.  Let  us  now  return  to  consideration  of  the  interaction 
correction  to  the  transmission  coefficient.  Since  this  is  calculated  from  the 
forward  far -field  amplitude  it  seems  that  both  kinds  of  simplification  we  have 
described  are  available.  The  transition  f rom  E quations  ( 3lio )  through  (3lb)  to 
(3hd)  can  now  be  explained  conceptually.  In  (3iic)  the  second  excitation  is  a 
line  source.  Its  strength  is  derived  from  interaction  consideration,  and  then 
slightly  simplified. 

In  (3iib)  this  line  source  has  been  effectively  replaced  by  the  equi- 
valent plane  wave .  The  influence  of  interaction  upon  the  strength  of  the 
secondary  excitation  is  still  retained.  The  same  is  true  in  (3hd)>  but  in 
addition  the  effect  of  interaction  on  the  strength  of  the  (simplified)  secondary 
excitation  has  been  obliterated;  the  influence  of  interaction  is  observed  only 
in  the  secondary  response.  Table  II  shows  the  wavelength  dependence  of  decrease 
in  accuracy  entailed  by  the  siirplifications  in  question.  Hence  from  this  table 
one  may  form  a  conception  of  the  relative  usefulness  of  the  various  idealizations 
we  have  been  discussing, 

V.   Excitation  by  a  line  so\irce 

When  the  slit  is  excited  by  a  line  source  instead  of  a  plane  wave,  an 
approxLmation  to  the  field  can  be  obtained  by  the  same  method  as  that  used  in 
Section  III,  except,  of  course,  that  now  each  half -plane  is  excited  by  two  line 
sources  instead  of  a  plane  wave  and  a  line  source.  We  append  here  a  brief  dis- 
cussion of  this  case. 

In  this  section  we  use  the  new  notation  E„  '(r  ,9  ,r,&)J.  for  the 
exact  response  field  produced  when  an  isolated  half -plane  is  excited  by  a  line 
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source  (r  ,9  ).  This  field  is  found  by  subtractinp;  from  the  total  field  (2)  the 
o  0 

corresponding  incident  field  (see  Fig, 2),  As  noted  earlier  ,  the  far  response 
field  can  be  expressed  asymptotically  and  uniformly  as  an  outgoing  wave  whose 
amplitude  is  a  function  of  (r  ,©,©).  The  (non-uniform)  asymptotic  expression 
for  the  diffracted  field,  on  the  other  hand,  is 

-ikr 

k 


(35)       [^l^^^^V^o'^'^^Ji"  V^  ^^;:r- ^(^o'®o*®)  ^°^  kr  »  i. 


Previously  we  have  used  the  function  f(r  ,6  ,9)  for  the  special  case    9  ■  n» 
Away  from  the  two  shadow  lines  9  »  n  *  9  this  function  is  given,  for  kr  »  1, 


by  the  asymptotic  form 

(35a)           f(r^,9^,9)^     -- 

1 

-i^o 
e 

■^  .  pr~. 

_ v  , 

COS  j(9-9q)   cos  5'(9+9^) 


as  can  be  shown  from  Equation  (2).  This  approximation  breaks  down  in  an  angular 

sector  surrounding  each  shadow  line  -  a  narrow  sector  if  kr  is  large  and  a  wide 

sector  if  kr  is  small, 
o 

Since  half -plane  (a)  is  excited  by  line  sources  at  (r  ,9  )  and  at 
(2d,n)  in  the  coordinate  system  associated  with  half -plane  (a),  while  half -plane 
(b)  is  excited  by  line  sources  at  (r,  ,9,  )  and  at  (2d,n)  in  the  coordinate 
System  associated  with  half -plane  (b),  the  total  field  E  is  given  by  the  sum 


(36) 


\''    p'"<'ao.«ao'--a'»a']i  *  [4"(--bo'V.--b'«b>]i  [  *  ^ 


xnc 
z 


E^'^^non)  +  F^'^^int). 
z  z 


*  Cf,  Eq.  (iia). 
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A  similar  equation  could  have  been  written  in  connection  with  Eq,  (13) 
of  Section  III.  Eq,  (36)  holds  approximately  for  any  point,  far  or  near,  as 
long  as  the  slit  is  wide  enouph.  We  can  now  evaluate  the  terms  of  (36)  for 
large  kr,  in  accordance  with  (35)  and  (35a),  and  get  an  equation  of  the  same 
type  as  Eq.  (13).  The  coefficients  C  and  C,  which  give  the  strength  of  the 

cL         O 

interaction  term  are  found,  by  the  method  of  Section  III,  to  be 

.         !^!W®bo»">  *  f(r^o»®ao»")^(2d,n,n) 
^a 


(37) 


1  -    (f(2d,n,n)]^ 
f(^o-.»8  .^'t)   +  f(r^^,a     ,n)f(2d,TT,n) 


-'  ^        ^   ao*  ao'   '             bo*  bo 
%~ 


1  -    |f(2d,n,7t)]' 


for  large  slit  widths  (kd  »  1), 

The  introduction  of  a  single  coordinate  system  such  as  employed  in 

Section  III  can  be  used  to  simplify  the  results.  A  further  simplification  ensues 

if  we  have  a  source  remote  from  the  slit,  i.e.,  if  kr  »  1.  In  addition  it  seems 

o 

useful  to  remark  that  the  calculation  of  far-fields  can  be  considerably  simplified 
by  the  use  of  the  reciprocity  theorem.  Even  in  the  presence  of  screens  the  latter 
theorem,  in  a  limiting  form,  can  be  used  to  express  the  far-field  of  a  line  source, 
in  terms  of  the  near  field  of  a  plane  wave.  In  consequence  of  our  approximations 
the  latter  field  for  the  slit  problem  can  be  expressed  in  terras  of  Fresnel 
integrals. 

Conclusion 

In  the  foregoing  work  we  have  derived  a  simple  and  convenient  expression 
for  the  field  diffracted  by  an  infinite  slit  in  a  conducting  plane,  i.e.,  when  the 
wavelength  is  greater  than  or  equal  to  the  slit  width.  The  principal  result  is 
that  this  field  can  be  accurately  calculated  everywlere  by  considering  each  half- 


-  32  - 


plane  composing  the  screen  to  be  excited  by  the  incident  field,  and,  in  addition, 
by  a  'virtual'  line  source  located  at  the  edge  of  the  opposite  half -plane  [Eqs. 
(13)  -(13e)3.  The  strength  of  this  line  source  can  be  simply  expressed  in  terms 
of  trigonometrical  functions  and  Fresnel  integrals  which  depend  upon  the  ratio 
of  slit  width  to  wavelength,  [Eqs,  (20),  (6)  and  (11)^.  These  expressions  may 
be  further  simplified  for  very  large  kd,  but  should  be  retained  intact  for  use 
with  slits  of  the  order  of  one  or  two  wavelengths  (cf.  conclusion  of  Section 
IV,  and  Table  if]. 

Up  to  the  present,  calculations  for  moderate  to  large  slit  widths  have 
required  the  use  of  the  Mathieu  fiinction  theory.  The  convergence  of  such  calcu- 
lations becomes  very  slow  as  the  slit  width  increases  to  three  wavelengths.  In 
addition  only  the  case  of  a  normally  incident  plane  wave  has  been  considered 
numerically  in  the  literature.  On  the  other  hand,  the  present  procedure  permits 
a)  a  rapid  calculation,  for  slit  widths  greater  than  one  wavelength j  b)  the 
accuracy  and  computational  convenience  increases  with  slit  width j  c)  the  pro- 
cedure is  capable  of  yielding  both  near  and  far  fields  in  terras  of  the  known 
solutions  for  a  half -plane j  d)  the  method  is  equally  well  applicable  to  more 
general  excitations  (e.g.,  oblique  plane  wave  or  line  source). 
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Table  I 


Legend:  The  first  column,  headed  F  (^),  gives  the  exact  angular  dependence 

of  the  far  field.   [Eq.  (32)1.  The  second  column,  headed  F   (0), 

'—  -i  '  non   ' 

gives  the  angular  dependence  calculated  from  the  non-interaction 
approximation.  The  third  column  gives  the  angular  dependence  calcu- 
lated from  the  approximation  developed  in  the  present  paper,  in  which 
interaction  between  half-planes  is  taken  into  account  [Eqs,  (3C),(3ir]< 
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Legend: 
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Table  II 
Transmission  Coefficients 

Exact  value  from  formula  (3Ua)  and  Table  I. 

Exact  value  from  Skavlem*-^^ 


T,    - 


Value  from  numerical  integration  of  exact  values  of  the  far  field, 
(cf.  first  column  of  Table  I) 

Value  from  numerical  integration  of  approximate  values  of  the  far 
field  (cf.  last  columr  of  Table  I) 


Tc;  ■  Value  from  formula  (3iic) 
Tz   ■  Value  from  formula     (3Ub) 
T»  =  Value  from  formula  (3Ud) 
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